Exact enumeration of the Critical States in the Oslo Model 
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We determine analytically the number Nn(L) of recurrent states in the Id Oslo model as a function 

of system size L. The solution N-jz(L) = + {^T^J is in exact agreement with 

the number enumerated in computer simulations for L = 1 — 10. For L>1, the number of allowed 
metastable states in the attractor increases exponentially as N-r(L) ~ c+A+, where A+ = 3+ 2^ is 
the golden mean. The system is non-ergodic in the sense that the states in the attractor are not 
equally probable. 

PACS numbers: 45.70.-n, 45.70.Vn, 05.65.+b 



Some non-equilibrium systems can be described by 
simple physical laws. Despite this being the case, we 
find complex and non-linear behavior in many such sys- 
tems. A slowly driven granular system illustrates this 
phenomenon spectacularly. If grains are deposited onto 
a finite surface, they will eventually build into a pile. The 
pile does not collapse and flatten due to the friction be- 
tween the grains. We call the different profiles of piles 
produced by this process, metastable states. The ad- 
dition of grains and the dissipation events (avalanches) 
form a delicate feedback loop. Avalanches initiated by 
dropping grains onto the pile will ensure that the slope 
does not become too steep, while the addition of grains 
prevents the slope from becoming too shallow. Even- 
tually, the system settles into an attractor where, the 
average avalanche size of the system diverges in response 
to small perturbations induced by adding grains. The 
distribution of avalanche sizes becomes scale invariant 
and a small perturbation may initiate small as well as 
large avalanches. As such, the system is intrinsically non- 
linear and displays characteristics similar to equilibrium 
systems poised at a critical point. This has been ob- 
served experimentally, for real granular systems where 
friction dominates over inertial effects 0, as well as nu- 
merically !J. It is an example of systems displaying 
self-organized criticality ]|, ||, ||. 

We study the attractor of a simple Id granular model 
known as the Oslo model. It successfully models slowly 
drivengranular systems displaying self-organized critical- 
ity m, |3) . The model has been shown to be a member of a 
large universality class of Id systems, which includes the 
de-pinning transition of an interface dragged through a 
random medium and the deterministic Burridgc-Knopoff 
train model for earthquakes H (7], g|. Thus, the Oslo 
model plays the role of the "Ising model" for Id self- 
organized critical s yst ems and has been widely studied 
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ever, in contrast to various deterministic models of granu- 
lar systems where exact results do exist ||| , the Oslo 
model has resisted analytic treatment || . Only a handful 
of exact results are known for non-equilibrium systems in 
general |23| and self-organized critical systems in partic- 
ular pi], B2l E4|. Any analytic result for the Oslo model 



would be very valuable. 

We are able to determine the number N-r (L) of states 
in the attractor of the_ Oslo Model. We show that 
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Nn{L) = i±# (3^1) — 
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sizes L » 1, the number increases exponentially fast as 
Nti(L) » c+A^, where A+ is the golden mean. This is a 
clear signature of the complexity displayed by this sim- 
ple model. Moreover, the system is non-ergodic. The 
recurrent states are not visited with the same frequency. 

The Oslo model is defined on a one-dimensional dis- 
crete lattice consisting of L sites, i = 1,2,... ,L. There is 
a vertical wall at the left boundary next to site i — 1 . The 
system is open at the right boundary next to site i = L, 
where grains are allowed to leave the system. We will de- 
note the respective heights of the pile starting from the 
vertical wall by an n-tuple h = (h\, h 2 , . . . , hi, . . . , hrj, 
see Fig. [l| In accordance with the experiment in Ref. [0 , 
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FIG. 1: The Oslo model of a Id granular pile. An integer 
height hi > is assigned to each site i. Grains are added 
at the site i — 1 next to the vertical wall by letting hi — > 
hi + 1. Grains can leave the system at the open boundary, 
hh — » hi, — 1. 

grains are only inserted at i = 1, the site next to the verti- 
cal wall. In the following discussion, it will be convenient 
to refer to local slopes instead of heights. The slopes are 
denoted by z = (z\ , z 2 , ■ ■ ■ , Zj, . . . , zr,) with the definition 
Zi = hi — hi + \ where /il+i = 0. Each site i is assigned a 
critical slope, z c (t) = (zf(t),z§(t),... ,zf (£),... ,z£(i)). 
If the local slope exceeds the critical slope, z, > zf(t), 
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a grain at site i will topple to the site i + 1, that is 
— * /ij — — * /ij+i + 1. As a consequence, the 
slopes at neighboring sites i±l may also exceed their crit- 
ical slopes and topple, consequently causing an avalanche 
that propagates until the system settles in a metastable 
state with Zi < zf for all i. 

The critical slope zf(t) is a function of time as it is 
chosen randomly every time site i topples. We denote 
the minimum allowed value of the random critical slopes 
z c(mm) anc j s j m ji ar iy the maximum allowed value z c ( max \ 
In general, we can use any discrete probability distribu- 
tion of critical slopes. However, we will consider the sim- 
plest case where zf is chosen at random between the 1 or 
2, zf(t) G {1, 2} when site i topples, such that z c ("»») = 1 
and z ci - max ^ = 2. 

The algorithm for the Oslo model is defined as follows. 

1. Initialize the critical slopes zf and the system in a 
arbitrary metastable state with Zi < zf for all i. 

2. Add a grain at site i = 1, i.e., z\ — > Z\ + 1. 

3. If Zi > zf, the site relaxes and 

Zi > Zi 2 
Zi±i -* z l± i + 1 

except for the sites at the vertical wall i = 1 and 
the open boundary i = L where 

Zi — ► Z\ - 2 z L — > z L - 1 

Z 2 — > Z2 + 1 — ► + 1- 

The critical slope at site i is chosen randomly if it 
topples, that is zf — > 1 or 2. A new metastable 
state is reached when < zf for all t. 

4. Proceed to 2 and reiterate. 

Note the separation of time scales which is built into 
the definition of the model. The addition of grains can 
only take place when the system has settled down into a 
metastable state. Thus the response of the system (the 
avalanches) is fast when compared with the interval be- 
tween perturbations. 

Metastable states M.j are rest states of the system. 
Different metastable states occur due to the different 
possible sets of critical slopes. Metastable states can 
also be divided into transient states Tj and attractor 
states IZj. Transient states are not reachable once the 
systems has entered into the attractor, which are recur- 
rent metastable states. The index j denotes discrete time 
steps associated with the long time scale of the system. 

When adding a grain to a metastable state Mj, it 
evolves into a new metastable state Mj+i by the relax- 
ation rules given above. Symbolically we write M.j > 
Adj+i, where the arrow is a shorthand notation for the 
operation of adding sand and, if necessary, relaxing the 
system until it reaches a new metastable state. Starting 



from, say the empty lattice T\ , which is a transient state 
as it will never be encountered again, we have 

l~i >—>■■■■ )—> T n i— > 7?4 i— ► ■ • ■ l— ► TZj-i i—* IZj i— > IZj+i 

After n additions and associated relaxations (if any), the 
system reaches the attractor of the dynamics. For exam- 
ple, if we set zf independent of time, we arrive at the 
BTW model in one dimension Any grain introduced 
into the system after the attractor state has been reached 
will simply flow to the open boundary fall out. In this 
case, there is only one attractor state with z, = zf Vi. 
We want to enumerate the number of different states in 
the attractor {IZj} for the Oslo model, where zf is chosen 
randomly between 1 or 2. 

We adopt the following nomenclature for sites 

Zl < z c( - mm ^> sinks, z t = z c{mm *> stable, 

z c( mm ) < z _ < z c( m ax) criticalj Zi > z c(max) supercriticaL 

Due to the way the problem is defined, Zj > always. 
Also, supercritical sites are by definition not allowed in 
any metastable state. Fig. |^ shows a diagrammatic re- 
lation of these classifications. From the relaxation rules, 

stable 

sinks sltes critical sites supercritical sites 

-) | ( 1 | ( 

z"""" 1 zf"" 1 

FIG. 2: We classify different types of sites according to the 
value of the local slope Zi. The above diagram shows these 
definitions in relation to one another, where z c i mm ) anc j 
z c(max) arg minimum and maximum allowed values for 
critical slopes respectively. 

we find that states in the attractor are subject to the 
following two constraints. 

(CI) Starting from the open boundary i = L, the first 
site (i < L) that is not a stable site (zi ^ z c(mm)^ 
must be a critical site. 

(C2) The first site that is not a stable state to the left 
of the sink must be a critical site. Equivalently, we 
cannot find a recurrent state in which two sinks are 
separated only by stable sites. The exception is at 
the closed boundary i = 1, where the first site that 
is not stable must be the boundary itself. 

We can show that these constraints must be fulfilled 
by assuming that each is not true and examining how 
that system would enter such as state. For example, Fig. 
||(a) shows a case where the the first rule CI does not 
hold. The possible preceding configurations are displayed 
in Fig. |^(b) and (c). Since these are not states in the 
attractor, the state displayed in Fig. ||(a) is not in the 
attractor. One could argue likewise for the second con- 
straint, C2. Note however, that by iterating these cases, 
we have only shown that these constraints are necessary 
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(a) (b) (c) 
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FIG. 3: (a) The configuration displayed violates the con- 
straint (CI), as the first site that is not stable (here, i = L— 2) 
next to the open boundary i = L is a sink, zl-2 — < 
2 c(mm) anc j ^ are p OSS j^] e preceding configurations. 
These configurations do not belong to the attractor. 



but not sufficient. However, we make the conjecture that 
these constraints are also sufficient for enumerating all 
attractor states from exact agreement with simulations 
for system sizes L = 1 to 10, see Fig. [|. 

For a system with L lattice sites, we want to find the 
number of states Nn(L) that obey the constraints above. 
We can simplify the problem further by first considering 
the number states with a fixed height at site i = 1, h\. 
We introduce the parameter q = hi — L = Ylf=i z i ~ L 
as it simplifies the equations involved considerably. Ob- 
serve that the range of q is < q < L. Let K(L, q) 
describe the number of states in a system of L sites 
and h\ = q + L. N-ji(L) can then be determined by 
Nn{L) =Y: L q=Q K{L, q ). 

To generate all states {z} that obey the counting con- 
straints for states in the attractor, consider the situation 
where we know all the allowed profiles for system sizes 
V < L. We denote the family of n-tuples that represent 
allowed profiles for a system of size L by J-(L). N-ji(L) 
is the number of n-tuples in ^(L). We can generate 
J-(L) — {h} by extending T(L — 1) by one site at the 
closed boundary. This may be done using the following 
procedure - extend each z' to z by writing Zi+i = z[. 
Consider one of three possible values for the new site 
next to the closed boundary, z\. zi = 2 for critical sites, 
zi = 1 for stable sites and zi = for sinks. Lastly, check 
that the new states in each of the three cases obey the 
counting constraints (CI) and (C2). 

Consider the three possible values for z%, that give 
J~(L) by extending T(L — 1) corresponding to three pos- 
sible cases such that K(L, q) — ^2 i=1 K^(L, q). 

1. When zi = 2, there is a critical site next to the 
closed boundary. The number of allowed profiles 
for a system of size L and parameter q must be 
the same as the number of allowed profiles with 
a system size of L' = L — 1. We relate q 1 to q by 
observing that q = hi — L — h[ — (L— 1)+1 = q' + l. 
Therefore we assert that when zi = 2, K^'(L, q) = 
K(L-l,q-l). 

2. Similarly, when zi = 1 corresponding to a stable 
site next to the closed boundary all profiles for a 
system size L — 1 with q' — q are allowed. In this 
case, (L, q) = K(L -l,q). 



3. Finally, the situation when z\ = means that 
a sink resides next to the closed boundary. The 
counting constraints state that the sink next to 
the closed boundary must be followed by a criti- 
cal site with any number of stable sites in between. 
In this case, we consider z' — (0, 2, ... ) extending 
T{L-2), z' = (0,1,2,...) extending T{L - 3) and 
so on to give L. This involves summing all allowed 
metastable states 0-(JJ) with L' < L~2. As above, 
we find that q' = q for all L'. Hence 

L-2 

K^(L,q) = Y / K( l ,q). 

i=l 

Note that this equation is a simple summation be- 
cause of the appropriate choice of q. 

Summing the three cases, 

L-2 

K(L,q) = K(L-l,q-l) + K(L-l,q) + J2K(i,q) 

i=\ 

L-l 

= ^K(i,q)+K(L-l,q-l) 

i=i 

= 2K(L -l,q) + K(L - 1, q - 1) - K{L -2,q-l). 

Hence we calculate the number of states as a function of 
system size by summing through q 

L 

N n (L) = J2 K ( L ^) 

9=0 

L-l L-2 

= 3j2K(L-l,q)-J2 K ( L -^l) 

q=0 q=0 

= W n (L-l)-N n (L-2). 

We have used the fact that < q < L. Nn(L) = 
3N-jz(L — 1) — N-jz(L — 2) is a Fibonacci like relation. 
Substituting a trial solution N-ji(L) = c+A^ + c_A^, 
and demanding the initial conditions Nn(l) = 2 and 
N-r{2) = 5, we obtain the solution, see Fig|j. 

N n{ L) = i±# ( 3 -±^ V + ^ V • 

K ' 2V5 \ % J 2V5 \ 2 J 

(1) 

The golden mean A+ = w 2.618 ... . Computer 

simulations were used to test these results for systems of 
size 1 to 10. They were found to be in full agreement. 

This result should be compared with the size of the 
attractor for the Abelian Bak-Tang-Wiesenfeld sandpile 
determined by Dhar p]] . He showed that N-ji = DetA, 
where A is the toppling matrix. In 2D and for large sys- 
tem sizes, N-ji increases exponentially fast with the num- 
ber of sites in the system (3.21 . . . ) LxL [p5| . Moreover, it 
was proven that the BTW model is erg odic |il). We find, 
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FIG. 4: Inset: Both simulation and analytic results are shown 
for systems with 1 to 10 lattice sites. Results obtained by 
enumerating the different metastable states in the attractor 
(filled circles) are identical to analytic results in Eq. (Q) (solid 
line). Main figure: The analytic solution of Eq. (|l|) showing 
the number of states in the attractor as a function of system 
size. For a moderate system size of L = 195, the number of 
metastable states in the attractor exceeds 10 81 , the estimated 
number of atoms in the universe. 

in contrast, that the Oslo model is non-ergodic in the 
sense that the states in the attractor are not visited with 
the same frequency. This is most easily demonstrated by 
showing that the probability of visiting stationary states 
in very small systems (e.g., L = 2 which has a total 



of 5 states in the attractor) is non-uniform. The num- 
ber of states in the attractor N-ji does not depend on 
the probability distribution with which z l c is chosen to 
be 1 or 2, except for the trivial cases of P(z l c = 1) = 
or 1, but the distribution will affect the accessibility of 
a given microstate. Finally, note that the Oslo model 
can be mapped into an interface dragged through a ran- 
dom medium ||. Thus, the results obtained for the Oslo 
model are directly applicable to the latter model as well. 

Possible extensions. A pile can be divided into an in- 
active zone consisting of grains that never will take part 
in the dynamics and the active zone. The width Xl of 
the active zone plays an important role as it is related 
to various physical phenomena. E.g. the average transit 
times of individual grains has been measured to increase 
with system size as (T) L oc LXl oc L 1+x . For the Oslo 
model x ~ 0-25 while experimentally \ ~ 0.5 @, |). We 
believe that using a similar approach one might be able 
to calculate analytically the distribution of the height of 
the pile at the vertical wall P{h\) and the critical ex- 
ponent X- This would be of great importance as this 
exponent will be the same in all the models belonging to 
the universality class of the Oslo Model. 
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